RESEARCH COMMUNICATION 



3D chromatin interactions 
organize Yan chromatin 
occupancy and repression 
at the even-skipped locus 

Jemma L. Webber, 1 Jie Zhang, 1 2 

Aaron Mitchell-Dick, 1 and Ilaria Rebay 1,2 ' 3 

^en May Department for Cancer Research, University of 
Chicago, Chicago, Illinois 60637, USA; 2 Committee on Cancer 
Biology, University of Chicago, Chicago, Illinois 60637, USA 

Long-range integration of transcriptional inputs is criti- 
cal for gene expression, yet the mechanisms remain 
poorly understood. We investigated the molecular de- 
terminants that confer fidelity to expression of the heart 
identity gene even-skipped (eve). Targeted deletion of 
regions bound by the repressor Yan defined two novel 
enhancers that contribute repressive inputs to stabilize 
tissue-specific output from a third enhancer. Deletion of 
any individual enhancer reduced Yan occupancy at the 
other elements, impacting eve expression, cell fate 
specification, and cardiac function. These long-range in- 
teractions may be stabilized by three-dimensional chro- 
matin contacts that we detected between the elements. 
Our work provides a new paradigm for chromatin-level 
integration of general repressive inputs with specific pat- 
terning information to achieve robust gene expression. 
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Historically, eukaryotic enhancers have been character- 
ized as DNA elements that are both necessary and 
sufficient to recapitulate endogenous gene expression. 
Many genes, particularly those with complex and dy- 
namic expression patterns, are controlled by multiple 
enhancers (Bulger and Groudine 2011), with the final 
pattern reflecting additive outputs from individual au- 
tonomously functioning elements. The advent of ge- 
nome-wide chromatin occupancy studies has revealed a 
far greater complexity to transcription factor-binding 
patterns than anticipated from this model, and the field 
remains in the early stages of deciphering the underlying 
functional relevance. For example, the contribution of 
nonadditive interactions between cis-regulatory elements 
to gene expression is just beginning to be appreciated, and 
the mechanistic importance of chromatin environment in 
organizing these interactions remains largely speculative 
(Frankel et al. 2010; Perry et al. 2010, 201 1 ; Dunipace et al. 
2011; He et al. 2011). 
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The transcriptional effector switch that controls cell 
fate specification downstream from Ras/MAPK signaling 
provides a tractable vantage point for exploring these 
ideas (O'Neill et al. 1994; Gabay et al. 1996; Graham et al. 
2010). The opposing action of two conserved ETS family 
transcription factors, the repressor Yan and the activator 
Pointed (Pnt), provides a bistable system that transitions 
from a high to a low Yan state to activate programs of 
gene expression that drive differentiation. For example, 
Yan and Pnt provide essential input to a combinatorial 
code that regulates expression of even-skipped (eve) during 
specification of the segmentally arrayed muscle and car- 
diac progenitors that arise in the dorsal mesoderm during 
embryonic stage 1 1 . While specification of these cells is 
reduced in pnt mutants and increased in yan mutants 
(Halfon et al. 2000; Knirr and Frasch 2001; Han et al. 2002), 
how the antagonistic transcriptional activities of Yan and 
Pnt are organized at the chromatin level is not understood. 

Here we tested the hypothesis that long-range integra- 
tion of Yan-mediated repression across a locus can stabi- 
lize the high Yan state to prevent inappropriate cell fate 
specification. Using regulation of mesodermal eve as the 
model, we found that targeted deletion of individual Yan- 
bound elements disrupts Yan occupancy across the locus 
and identified a three-dimensional (3D) chromatin con- 
formation that may stabilize these long-range interactions. 
Loss of coordinated Yan-mediated repression compromises 
the Yan-Pnt switch, as evidenced by increased variability 
of eve expression and impaired cell fate specification, 
cardiac function, and viability — all phenotypes that can 
be enhanced or suppressed by reducing the genetic dose of 
yan or pnt, respectively. We speculate that chromatin- 
level integration of repressive inputs with specific pattern- 
ing information may provide a broadly used mechanism to 
achieve robust gene expression. 

Results and Discussion 

In addition to the muscle/heart enhancer targeted by the 
Yan-Pnt switch (MHE) (Halfon et al. 2000), our genome- 
wide occupancy study of Yan mapped three other regions 
to the eve locus (Webber et al. 2013). We refer to the four 
Yan-bound domains as Dl, D2, D3/MHE, and D4. To 
explore the contribution of individual Yan-bound regions 
to mesodermal eve expression we recombineered specific 
deletions into a functional genomic bacterial artificial 
chromosome (BAC) construct encompassing the entire 
16.4-kb eve locus carrying an in-frame YFP tag (Fig. 1A ; 
Ludwig et al. 2011). As enhancer-blocking activity has 
previously been identified between the MHE and Ter94 
gene (Fujioka et al. 2009), we reasoned that the D4 Yan- 
bound element was more likely to contribute to Ter94 
than to eve regulation and excluded it from our analysis. 
Transgenes were generated, and altered function was 
assessed by examining expression in the Eve-positive 
cardiogenic precursors. The simple prediction for deletion 
of a functional Yan-bound region is that loss of repression 
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Figure 1. Yan-bound regions at the eve locus contribute to meso- 
dermal expression, viability, and cardiac function. (A) Yan ChlP-chip 
patterns at eve. Yan binds to the proximal promoter (Dl), a region 
that partially overlaps a cis-regulatory module responsible for 
neuronal and anal plate expression (D2), the previously character- 
ized MHE (Half on et al. 2000) (D3), and a region upstream of the 
TER94 locus (D4) that overlaps a previously characterized insulator 
and homing element (Fujioka et al. 2009). Red boxes depict engi- 
neered deletions. Gray boxes depict stripe enhancers. Light-orange 
boxes depict regions used for reporter constructs. (B-E) Immuno- 
staining of eve WT (B), eve AD1 (C), eve Aul (D), and eve ADS (£) with anti- 
GFR Eve-YFP expression recapitulates wild-type eve expression in 
1 1 clusters of cells in the mesoderm (blue arrowhead) and the CNS 
(red arrowhead). While eve Ar>1 and eve Ar>1 qualitatively reproduced 
eve WT expression, reduced expression was observed with eve Ar>3 . 
(F) eve AD , eve Ar>1 , and eve A embryos were stained and fixed in 
parallel with eve WT embryos. Quantification of pixel intensity for 
each cluster relative to eve WT is depicted on the Y-axis. Mean 
intensity is indicated with a red line, and error bars indicate standard 
deviation, (n > 30 embryos; eve WT :eve Ani , P < 0.0001; eve WT :eve AJD1 , 
P < 0.0001; Mann-Whitney U ; n = 12 embryos; eve WT :eve ADS , 
P < 0.0001). (G) eve-nuU embryos [Di(2R)eve/eve R13 ] carrying the 
deletion transgenes were genotyped and cultured at 25 °C, and the 
number of adults that survived was counted [n > 350; eve WT :eve AJD1 , 
P < 0.0001; eve WT :eve AD1 , P < 0.0001; eve WT :eve ADS , P < 0.0001). 
(H) Average heart rates of white prepupae. (eve WT :eve AD1 , P < 0.0001; 
eve WT :eve Ar>3 , P < 0.0001; Mann- Whitney U). (I) In the flight assay, 
a flight index of 4 is the highest possible score; 0 indicates no 
flying ability (eve NT \eveF Dl , P < 0.01; eve^-.eve^ ', P < 0.0001; 
Mann- Whitney U). (/) eve Ar>1 and eve Ar>1 animals had a reduced 
negative geotaxis response compared with eve WT (n > 100). 



should increase expression of the target gene,- however, if 
the deletion also removes binding sites for critical activa- 
tors, then reduced expression is expected. Either scenario 
predicts reduced ability to complement an 6V£-null. Alter- 
natively, if elements are functional yet fully redundant 
under optimal conditions, no change in expression or 
viability is predicted. 

The control eve-YFP transgene fully complemented an 
£ve-null background, with 86% rescue to adulthood, and 
was expressed in a pattern identical to that of endogenous 
eve. Thus, in stage 1 1 embryos, Eve-YFP was prominent 
in the muscle and cardiogenic precursors as well as in 
clusters of neurons in the CNS (Fig. IB). Deletion of the 
Dl, D2, or D3 Yan-bound regions altered expression and 
reduced genetic rescue efficiency (Fig. 1C-G), indicating 
important contributions to regulation of eve. As predi- 
cted based on loss of key activating inputs (Halfon et al. 
2000), deletion of the D3/MHE region significantly re- 
duced mesodermal expression, with Eve AD3 -YFP expres- 
sion detectable, on average, in only one cell per cluster 
(Figs. IE, 2A). For the Eve AD1 -YFP and Eve AD2 -YFP trans- 
genes, although expression in the mesoderm appeared 
qualitatively normal (Fig. 1B-D), quantification of the 
YFP signal revealed significant increases in both the 
mean intensity and embryo -to -embryo variation relative 
to the wild- type Eve-YFP control (Fig. IF). The reduced 
ability of the deletion transgenes to rescue an 6V6-null 
(Fig. 1G) suggests that these Yan-bound elements are not 
redundant. The deleted regions do not include known 
stripe enhancers (Fig. 1A, gray boxes), and cuticle prepa- 
rations of dead embryos from the rescue experiments did 
not show segmentation defects (data not shown). Thus, 
the reduced viability is unlikely attributed to axial pat- 
terning defects. 

Precise regulation of eve expression in the cardiogenic 
mesoderm is essential for acquisition of heart cell iden- 
tity and ultimately for cardiac function (Fujioka et al. 
2005; Johnson et al. 2007). Consistent with this, the 
elevated expression of Eve AD1 -YFP reflected in part an 
increase in the number of Eve-positive mesodermal cells 
specified (Fig. 2A), raising the possibility that compro- 
mised heart development might contribute to reduced 
fitness. To test this, we first measured pupal heart rates 
(Fig. 1H). In agreement with Fujioka et al. (2005), the heart 
rate of eve~'~; £V£ AD3 -YFP animals was significantly 
reduced relative to the eve~^~; eve-YYV control. eve~^~; 
£ve AD1 -YFP larvae also had a reduced heart rate, suggest- 
ing that the ectopic cell fate specification detected 
compromises the ensuing developmental program. De- 
letion of the D2 region did not significantly change the 
number of Eve-positive mesodermal cells specified, and 
the heart rate measured in eve~ ! ~ 
parable with that of the control. 

Adult heart function is commonly assessed with assays 
that measure the general vigor and activity of the fly 
(Fujioka et al. 2005; Piazza et al. 2009; Shahrestani et al. 
2009). In both a flight assay and an assay measuring the 
negative geotaxis response before and^ after heat stress, 
eve~!~; ev6 AD3 -YFP and eve~ l_ ; eve 
performed poorly relative to eve~^ 
animals (Fig. II J). eve~^~; 6V6 AD2 -YFP animals exhibited 
a response intermediate to that of the eve~^~; 6V£ AD1 -YFP 
and control eve~^~; eve-YVV animals in the geotaxis assay 
prior to heat stress. However, after stress, the activity of 
eve~/~; £V£ AD2 -YFP flies was significantly below that of 



eve Ar>2 pupae was com- 



AD1 -YFP flies again 
eve-YW control 



2294 GENES & DEVELOPMENT 



Long-range Yan-mediated repression of even-skipped 




Eve WT Eve AD1 Eve AD2 Eve AD3 



4* ^ ^ f f & 

,^ .4® /A® ,4® 

<c/ .V V <<? .<y <<? .<y 




C CO ' 



eve-/-; Eve WT 
• eve-/-; Eve AD1 
eve-/-; Eve AD2 



eve-/-; Eve WT, pntAF397/ 
+, Eve WT 

eve-/-; Eve AD1, pntAF397/ 
Eve AD1, + 

eve-/-; Eve AD2, pntAF397/ 
Eve AD2, + 



123456789 41424344454647484950 

Time (min) 

Figure 2. Robust regulation of mesodermal Eve expression depends 
on yan and pnt. [A) Eve-positive mesodermal cells were counted in 
eve WT , eve Ar>1 , eve Ar>1 , and eye AD3 embryos in the presence or 
absence of t wis t-Gal4>U AS Yan ACT . A statistically significant in- 
crease was observed in eve AID1 , and a reduced number was observed 
in eve AD3 (n > 70 embryos; P < 0.0001 and P < 0.0001, respectively; 
Mann- Whitney U). [B] Quantification of Eve-YFP shows an increase 
in intensity and variation of expression in embryos heterozygous for 
yan (n > 10 embryos; P < 0.0001; ANOVA). Heterozygosity for yan 
further increased YFP intensity in eve*' 02 but not eve Ar> [n > 10 
embryos; P < 0.0001; ANOVA). (C) Hand-genotyped eye-null em- 
bryos [Di{2R)eve/ eve R13 ] or eye-null, yan heterozygous embryos 
[Df(2R)eve, yan 833 /eve ] carrying the different deletion transgenes 
were cultured at 25°C, and the number of eclosed adults was counted 
(n > 250; eve WT :yan 833 /+ ; eve WT , P = 0.0007; eve AD3 :yan 833 /+; 
eve Ar>3 , P = 0.0013). (D) Quantification revealed an increase in 
Eve-YFP intensity in eve Ar>1 and eve Am embryos relative to eve WT 
(P < 0.0001). No significant difference was observed between eve AE>1 
or eve AID2 and eve WT Eve-YFP expression in pnt heterozygotes. [E) 
Flies of the indicated genotype were subjected to a geotaxis re- 
sponse assay. Heterozygosity for pnt restored the response of 
eve AD1 to that of an eye WT control [n > 50) and improved that of 
e^e Am (n > 35). 



the control. Together, these results confirm the physio- 
logical importance of the Yan-bound Dl region and 
support the hypothesis that it mediates repressive inputs 
important for stabilizing mesodermal eve expression 
within levels required to support normal cell fate acqui- 
sition during heart development. As the D2 deletion 
partially overlaps a neuronal enhancer, further analysis 
will be required to determine whether the reduced fitness 



and activity result from either compromised neural 
function, subtle cardiac defects, or a combination of the 
two. 

To obtain more direct evidence that compromised Yan- 
mediated repression contributes to increased expression 
and reduced fitness of the deletion mutants, we first 
tested the effectiveness of twist-Gal4- driven overexpres- 
sion of a constitutively active form of Yan, Yan ACT (Rebay 
and Rubin 1995), in repressing expression of the deletion 
transgenes. If Yan-mediated repressive inputs act exclu- 
sively on the MHE, then Yan ACT should repress mesodermal 
expression of Eve AD1 -YFP and Eve AD2 -YFP but not Eve AD3 - 
YFP. Alternatively, if Yan-mediated regulation across 
multiple genomic elements contributes to stable eve 
expression, then all deletion transgenes should be sensi- 
tive to Yan ACT . Our results support the latter scenario, as 
Yan ACT fully repressed all three eve-YW deletion trans- 
genes (Fig. 2A). 

We next asked whether heterozygosity for yan or pnt 
might, respectively, enhance or suppress the loss of 
robustness measured for the eve deletion constructs. 
Using either a null allele or a small deficiency, we found 
that heterozygosity for yan increased the mean intensity 
and variation of wild-type Eve-YFP expression and further 
increased the intensity of Eve AD2 -YFP expression (Fig. 
2B). While yan heterozygosity did not significantly in- 
crease Eve AD1 -YFP expression, the number of Eve-posi- 
tive cells specified was increased (Supplemental Fig. SI). 
These changes also impact fitness, as yan heterozygosity 
further reduced the genetic rescue efficiency of the eve 
transgenes (Fig. 2C). Conversely, heterozygosity for pnt 
suppressed the increased intensity and variability of 
Eve AD1 -YFP and Eve AD2 -YFP expression (Fig. 2D ; Supple- 
mental Fig. SI). A reduced pnt dose also improved the 
negative geotaxis response following heat shock (Fig. 2E). 
Together, these results suggest that the Dl and D2 
elements are required for robust regulation of eve expres- 
sion in the muscle and cardiogenic precursors such that 
loss of Yan-mediated repression upon deletion of either 
element permits inappropriate Pnt-mediated activation 
of eve, presumably through the MHE, leading to in- 
creased and more variable expression. Such destabiliza- 
tion of the Yan-Pnt switch makes the system prone to 
aberrant cell fate specification and ultimately reduces 
fitness. 

To explore further the properties of the Dl and D2 Yan- 
bound regions, we asked whether they define autono- 
mous cis -regulatory elements or whether their activity 
requires the endogenous genomic context. Reporter con- 
structs carrying sequences spanning the chromatin im- 
munoprecipitation (ChlP)-defined Yan-bound regions 
were generated (Fig. 1A), and expression was examined 
in stage 11 embryos. In contrast to the control MHE 
reporter that expressed specifically in the expected me- 
sodermal clusters, the Dl and D2 reporters showed weak 
ubiquitous expression (Fig. 3A). A slightly stronger signal 
was detected in the mesoderm for both reporters but in 
domains broader than the three cell clusters seen with the 
MHE reporter. Removing yan increased the intensity but 
not the specificity of the expression patterns of the Dl 
and D2 reporters, suggesting that the individual elements 
can be recognized and repressed by Yan outside of the 
context of the intact eve locus. The autonomous ability of 
each element to recruit Yan was confirmed by ChlP-PCR, 
although quantitative PCR (qPCR) analysis revealed that 
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Figure 3. Regions Dl and D2 are sufficient to recruit Yan but not to 
drive specific mesodermal expression. [A) Analysis of GFP expres- 
sion of the Dl, D2, and D3 reporters (depicted in Fig. 1A) in a wild- 
type or yarT^ (yan ER433 ) background; staining and imaging of each 
pair was done in parallel. The D3/MHE reporter expressed strongly 
and specifically in the mesodermal clusters. The Dl and D2 reporter 
expression was weak and ubiquitous with only slight mesodermal 
enrichment. Expression intensity increased for all reporters in yan~^. 
(B) ChlP-PCR shows Yan enrichment at the Dl, D2, and D3 
reporters. To the distinguish reporter from the endogenous locus, 
primer pairs included one within the reporter and one from the 86Fb 
integration site. Enrichment at argos (aos) and mae provided positive 
controls. (C) ChlP-qPCR showed less efficient Yan recruitment to 
the reporters relative to the endogenous eve locus. Data are plotted 
relative to Yan occupancy at aos. Error bars are mean ± standard 
deviation (n = 4). 



occupancy was reduced relative to that measured at the 
endogenous locus (Fig. 3B,C). 

Because the Dl and D2 Yan-bound elements do not 
appear to provide independent additive inputs capable of 
directing mesodermal eve expression, we wondered 
whether their primary role might be to integrate Yan- 
repressive inputs across the intact locus. We hypothe- 
sized that long-range interactions between these regions 
might stabilize Yan binding, thereby providing a buffering 
mechanism that sets the precise MHE-driven eve expres- 
sion pattern needed for accurate cell fate specification. If 
correct, then deletion of any individual Yan-bound region 
not only should abolish occupancy at the site of deletion, 
but might also perturb binding at the remaining elements. 
To test this, we assessed Yan occupancy by ChlP-qPCR in 
embryos homozygous for the eve deficiency and carrying 
the deletion transgenes. Supporting the model of coordi- 
nated stabilization of Yan occupancy, deletion of any 
individual Yan-bound domain reduced Yan occupancy by 
about twofold at the remaining two eve elements but did 
not alter enrichment at a distant locus (Fig. 4A,B; Supple- 
mental Fig. S2). Thus, precise regulation of mesodermal 
eve expression relies on long-range coordination between 
cis-regulatory elements rather than additive inputs from 
autonomous enhancers. 

One possible mechanism is that the 3D chromatin 
environment facilitates interactions between the Dl, D2, 
and D3 elements. Using proximity ligation and PCR 
analysis (referred to as chromatin conformation capture 
[3C]), we detected a long-range interaction between a Dl- 
D2 fragment and the D3/MHE fragment (Fig. 4A,C,D). 
Sequencing confirmed the expected identity of the am- 
plified ligated fragment. Control primers located just 
outside of the eve locus did not yield a 3C product with 
the D1-D2 fragment (Fig. 4D ; Supplemental Fig. S3A,B). 



Furthermore, the 3C interaction was not detected in 
adult heads, suggesting context-specific chromatin con- 
tacts (Supplemental Fig. S3C-E). Finer-scale mapping in 
both wild-type and yan mutant embryos should provide 
further insight into how the 3D chromatin environment 
facilitates cooperative recruitment of Yan across the 
locus and how this in turn both stabilizes the Yan-Pnt 
switch in optimal conditions and buffers against genetic 
or environmental variation that might limit Yan. 

In conclusion, our results suggest a novel regulatory 
mechanism in which, rather than contributing patterning 
information, the Dl and D2 enhancers act to dampen 
the expression control mediated by the D3/MHE. More 
broadly, this raises the possibility that even when tran- 
scription factor-bound regions identified through ChIP 
sequencing (ChlP-seq) fail to give clear expression pat- 
terns as autonomous reporters, they may still contribute 
to gene regulation in vivo. In the case of Yan, we speculate 
that its intrinsic polymerization ability (Qiao et al. 2004; 
Song et al. 2005; Zhang et al. 2010) allows it to exploit the 
3D chromatin environment to organize long-range re- 
pressive complexes that coordinate information across 
multiple enhancers. 
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Figure 4. Yan binds cooperatively to individual regions across the 
eve locus. (A) Depiction of the Yan ChIP pattern at the eve locus. 
Red boxes indicate regions deleted, blue dashed lines depict EcoRl 
sites, and arrows indicate PCR primers used for 3C analysis. (B) 
ChlP-qPCR shows enrichment was lost at the site of deletion 
(arrow). A consistent decrease in enrichment at the other Yan-bound 
regions was observed in all deletions relative to eve WT . Error bars are 
mean ± standard deviation. (C) Forward and reverse primers located 
in the Dl region provided a PCR-positive control (pc), amplifying 
a fragment in both cross-linked and uncross-linked DNA. A PCR 
product from a primer located at the end of the fragment containing 
the Dl and part of the D2 region (primer 1; light blue arrow) and 
a primer located at the 3' end of the MHE (primer 2) fragment was 
observed only in the cross-linked sample. Asterisk indicates primer- 
dimer. A PCR product from primer 1 and a primer located at the 3' 
end of a fragment containing the neuronal enhancer but not the 
insulator/homing element (primer 3) was also observed in the cross- 
linked sample. Sequencing confirmed the predicted ligated products. 
(D) ChlP-qPCR of cross-linked and uncross-linked samples. The 
relative amount of 3C product was determined for primer pairs 1+2, 
1+3, and 1 + nc (nc primer located to the right of the insulator 
element) (see Supplemental Fig. S3A), and values were plotted 
relative to a positive control. Error bars are mean ± standard 
deviation (n = 3; cross-linked 1 + 2:uncross-linked 1 + 2, P = 
0.0015; cross-linked 1 + 3:uncross-linked 1 + 3, P = 0.0007; cross- 
linked 1 + nc:uncross-linked 1 + nc, P = 0.2723). 
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Many parallels can be drawn between Yan/Pnt regula- 
tion of the transition from uncommitted progenitor to 
specified cell fate within a developing tissue and the 
control of the pluripotent state in mammalian cells. 
Given the potential for natural, genetic, or environmental 
variation to perturb the timing, accuracy, or stability of 
such transitions, mechanisms to impart robustness to 
gene expression are almost certainly required and very 
likely conserved across species. Thus, the regulatory para- 
digm established in this study may be broadly relevant. 

Materials and methods 

Deletion transgenes and GFP reporters 

A BAC containing the eve locus (Ludwig et al. 2011) was modified to 
generate the Dl, D2, and D3 deletions using the Counter-Selection BAC 
Modification kit (Gene Bridges) following the manufacturer's protocol. 
The genomic coordinates for the Eve deletions Dl, D2, and D3 were 2R: 
5866200-5866500; 2R: 5869000-5869800; and 2R: 5872789-5873131, 
respectively. Subcloning via the BAC Subcloning kit (Gene Bridges) into 
the transformation vector and site-specific integration in the attP2 site on 
the third chromosome (Groth et al. 2004) were performed as described in 
Ludwig et al. ( 201 1). Reporter constructs were generated by inserting the 
Dl (2R: 5865760-5867730), D2 (2R: 5868845-5870866), and MHE (2R: 
5872818-5873045) sequences into the pJR20 plasmid upstream of a min- 
imal promoter and the GFP gene coding sequence. Site-specific integra- 
tion into the 86Fb site on the third chromosome was as described in 
Bischof et al. ( 2007). 



quadrants, with the highest denoted as 4 and the lowest denoted as 1 . Flies 
that fell through without touching the sides were scored as 0. 

Negative geotaxis response 

Assays were performed essentially as described in Tong et al. ( 2007) and 
Shahrestani et al. (2009). Ten 2- to 3-d-old males per genotype were 
collected into individual culture vials, and a horizontal line was drawn 
1 cm above the food. Flies were knocked to the bottom, and after 1 5 sec, the 
number that climbed above the 1 -cm line scored. This was repeated every 
30 sec for 10 min to establish a baseline response. Animals were then 
subjected to a 30-min heat shock at 37°C and returned to room tempera- 
ture, and the geotaxis response was determined as they recovered. This 
assay was repeated until n-values of >100 were achieved for eve-YFP, 
eve AD1 , eve AD2 , and eve AD3 ; >50 were achieved for eve-YFP, pnt, and eve AD1 , 
pnt } and >35 were achieved for eve AD2 and pnt. 

ChIP 

ChlPs from stage 11 embryos (5h20-7h20 at 25°C) were performed as 
described in Webber et al. ( 2013). Analysis of Yan binding to the Dl, D2, 
and MHE GFP reporters was performed with a primer specific to the pJR20 
plasmid and a primer specific to the Dl, D2, or MHE sequence. ChlP- 
qPCR was performed using the QuantiTech SYBR Green PCR kit 
(Qiagen). A standard curve for each primer pair was generated using serial 
dilutions of genomic DNA. The relative amounts of input, mock-treated, 
and immunoprecipitated DNA were determined based on the standard 
curves, and the ChIP signals were calculated as immunoprecipitation/ 
mock ratios. Each individual mock-treated and immunoprecipitation 
sample was derived from the same set of cross-linked embryos. 



Fly strains and genetics 

eve~!~; eve-YFP, eve _/ ~; eve AD1 , eve~ 7 ~; eve AD2 , or eve~'~; eve AD3 embryos 
were selected from the cross of Df(2R)eve (Bloomington Stock Center, no. 



eve R13 /CTG ; eve\ 



WT AD1 AD2 



or AD3 ). CTG-negative embryos heterozygous 



for yan and homozygous for the eve deletion of interest were selected from 
the cross between either yan 833 /CTG ; eve( WT , AD1 , or AD2 ) or Di(2L)yan/ 
CTG ; eve( WT or AD1 ) and +/+, eve( WT , AD1 , or Am ). 

To analyze eve-YFP expression in embryos heterozygous forpzzt, Ubx-lacZ- 
negative embryos were selected from a cross of pnt^ 397 , eve^, AD1 , or AD2 )/ 



TM3, Ubx-lacZ and +/+; evey 



Alternatively, TTG- 



negative embryos were examined from a cross of pi2£ A88 /TM3, twist- 
Gal4, UAS-GFP (TTG), and eve WT , eve AD1 , or +/+. To determine the 
effect of yan and pnt heterozygosity on the ability of eve[ WT , AD1 , or AD2 ) 
to rescue an eve-null, CTG-negative embryos were hand-selected from 
a cross of Df(2R)eve, yan 833 /CTG ; eve( WT , AD1 , AD2 , or AD3 ), and eve R13 / 



CTG; eve\ 



WT AD1 AD2 



or J or Df(2R)eve/CTG; pnt Am9J , evef 



and eve R13 /CTG; evef 



For ChlP-qPCR at the eve 



locus, -400 CTG-negative embryos were selected from the progeny of 



3C 

Cross-linking was performed as described in Webber et al. ( 2013) but for 
20 min. Embryos were homogenized in buffer A (15 mM HEPES at pH 7.6, 
10 mM KC1, 5 mM MgCl 2 , 0.1 mM EDTA, 0.5 mM EGTA, 350 mM 
sucrose, 1 mM DTT) and centrifuged for 15 min at 10,000 rpm. The pellet 
was resuspended in 1.5X EcoRI buffer (New England Biolabs) with BSA 
and incubated for 30 min at 37°C with 3 |xL of 10% SDS. Following 20 min 
at 65 °C, 20 |jlL of 10% TX-100 was added, and samples were incubated 
for 1 h at 37°C. Four microliters of EcoRI (20,000 U/mL ; New England 
Biolabs) was added and incubated overnight at 37°C and digestion was 
stopped by adding 16 |jlL of 10% SDS and incubating for 20 min at 65°C. 
Digested chromatin was diluted in ligation buffer (50 mM Tris at pH 7.8, 
20 mM DTT, 1 mM ATP, 25 |xg/mL BSA) with 1% TX100 to quench the 
SDS. Following 1 h of incubation at 37°C, 2 |xL of T4 DNA ligase (400,000 
U/mL; New England Biolabs) was added, and samples were incubated 
overnight at 16°C. Three microliters of Proteinase K (Invitrogen) was 
added prior to incubation overnight at 65°C. DNA was purified by phenol/ 
chloroform extraction and ethanol precipitation, and ligation products 
were analyzed by PCR or qPCR. 



Quantification of Eve-YFP expression 

Stage 11 embryos were stained as described in Webber et al. ( 2013) with 
rabbit anti-GFP (1:1000; Invitrogen). Using a Zeiss 510 confocal micro- 
scope, serial 0.8-|xm z-sections were taken through the Eve-positive 
mesodermal cells, and maximum was projected using LSM software. 
Expression intensity was calculated as the mean pixel intensity for the 
cluster minus the mean background pixel intensity. Measurements were 
normalized to eve WT . Column statistics are shown in Supplemental Table 1 . 

Heart rate measurement 

The number of heartbeats of white prepupae in a 15 -sec period was determined 
as described in Cooper et al. 2009) and charted as hertz (cycles per second). 

Flight assay 

Flies were tapped into a 1000-mL graduated cylinder coated with mineral 
oil as described in Cauchi et al. (2008). The cylinder was divided into four 



Statistical analyses 

Analyses were performed with Prism 5. A D'Agostino-Pearson normality 
test was used to determine whether the data fit a Gaussian distribution. A 
t-test or an ANOVA was performed on normally distributed data. 
Alternatively, a Mann- Whitney U or a Kruskal-Wallis test was used. A 
Dunn-corrected P-value is reported for statistical tests carried out on 
multiple data sets. The a level for all tests was 0.05. 
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